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It is shown that a superposition of static and rapidly oscillating electric monopole (source) fields 
is capable of trapping particles with a permanent electric dipole moment. Thus, the new trapping 
mechanism differs fundamentally from saddle-point traps that use static and oscillating higher- 
multipole fields. An analytical stability analysis together with detailed molecular dynamics trajec- 
tory calculations prove that the trap is stable. Thin rods of barium titanate (BaTiOs) provide an 
illustrative example for the working principle of the new trap. The effects of gravity are considered. 
. The existence of a bifurcation regime is predicted. A particular strength of the new trap is that it 

^S) ' also works for zero orbital angular momentum with respect to the field-generating electrodes. 



PACS numbers: 37.10.Pq, 77.22.-d 
! I. INTRODUCTION 



Levitation of microscopic and macroscopic objects is of general interest in physics: It is a key technology for basic 
^ and applied sciences. For instance, research areas ranging from high-resolution spectroscopy [H-Qi frequency standards 
'p^. [3| , neutral anti- matter production Q , Bose- Einstein condensation 0, |1] , and quantum computing [§| to high-speed 
trains [l^, all rely on the levitation and three-dimensional confinement of objects by electromagnetic fields. The key 
point is that electromagnetic confinement provides "walls of pure energy" that avoid contamination, friction, and, in 
the case of antimatter, annihilation, by avoiding contact with material walls. 
^ ' Traps used in science and technology may broadly be divided into charged-particle traps and neutral-particle traps. 
. Historically, charged-particle traps were developed first, starting with the development of the Kingdon trap [ll| in 
the 1920s, the Penning trap [l3 in the 1930s and the Paul trap 13] in the 1950s. In addition to their use in high- 



> 
1^ 



^ . resolution spectroscopy (H-y], freq uency standards Q, and quantum computing £], these traps were used, e.g., as a 



Oh' 



model for planetary dynamics [lj| , for studying collisions and excited-state lifetimes of ions [15| , and the production 
and investigation of Coulomb crystals [16|-[l8j. Recently, neutral-particle traps were developed and demonstrated as 
well. Based on the type of field ernployed, we may distinguish several different types of neutral-particle traps, for 
instance, electric [l9l - [2a | , magnetic [2^ 113] , or light-field traps [H, 111] . 

The purpose of this paper is to introduce a conceptually new trap for the stable confinement of neutral, polar 
particles. The trapping mechanism of the new trap is fundamentally different from the trapping mechanisms employed 
. in existing neutral-particle traps: (i) The trap does not make use of energy shifts in selected quantum states and (ii) 
' in contrast to existing traps |23l . |24| , which are based on a dynamic saddle-point potential akin to the strong- focusing 
^NJ ; Paul-trap principle [l3|, the new trap works on the basis of a superposition of static and oscillating monopole fields. 

The paper is organized in the following way. In Sec. |ll] we derive the coupled radial and angular equations of motion 
, of point dipoles trapped in a combination of ac and dc electric monopole (source) fields. In Sec. IIIII we analyze the 
I ■ stability of the coupled system of equations and derive conditions for the stable confinement of point dipoles. We 
conclude that stable trapping is possible as long as stable solutions of the radial equation exist. Therefore, in Sec. IIVI 
we investigate in detail the stability properties of the radial equation. We find that subject to certain conditions 
on the control parameter of the radial equation, the radial equation is stable, which guarantees the stability of the 
combined radial and angular system of equations derived in Sec. |IT1 and thus results in stable confinement of point 
dipoles according to the new trapping scheme. In Sec. |Vl strengthening our analytical results, we show numerically 
that the trap is stable. In Sec. IVII we derive the conditions under which stable confinement in the presence of gravity 
is achieved. In Sec. IVIII we state a step-by-step algorithm according to which trap parameters may be designed. The 
trapping of thin, cylindrical rods of barium titanate is discussed in Sec. I Villi In Sec. IIXI we discuss our results. In 
Sec.|X]we summarize our results and conclude the paper. 



II. EQUATIONS OF MOTION 



In this section we derive the equations of motion of a point dipole p of mass M and moment of inertia / in an 
electric field 

Eir,t) . ^r, (1) 



2 



where u > and a > are constants, r is the distance of the dipole from the source of the field, r is the unit vector in 
r direction, and V{t) is the voltage applied to the electrodes that generate the field E. Since the interaction energy 
of a point dipole p with an electric field E is 

W ^ -p E = -pEcosie), (2) 

where 9 is the angle between p and E, the Lagrangian function C of the point dipole in the case of zero orbital angular 
momentum is 

£ = hlr' + ll[e' + sin'i0)^'] + cos(0), (3) 

where ip is the azimuthal angle of p. Since dC/d(p — 0, the azimuthal angle is a cyclic variable [30l|, and the conjugate 
canonical momentum, 

- g - I sin' {9)^, (4) 

is a constant of the motion. In the following we specialize to the case L^p = 0. While some neutral-particle traps 
require non-zero orbital angular momentum to work properly (an example is the trap described in (l^jUll), a particular 
strength of the trap introduced in this paper is that it works for zero orbital angular momentum. Therefore, playing 
to the strength of the new trap, we treat the case in which the dipole is not rotating around the field source, i.e. the 
orbital angular momentum of the point dipole around the field-generating electrodes is zero. It has been checked by 
means of extensive numerical simulations that the trap also works for non-zero orbital angular momentum. However, 
since this case does not add substantial additional insight into the fundamental working principle of the new trap, the 
case of non-zero orbital angular momentum is not treated in this paper and we will continue to focus on the case of 
zero orbital angular momentum as reflected in the Lagrangian function ([3]) . On the basis of the La gran gian function 
^ the radial (r) and angular (9) equations of motion are obtained from the Lagrangian equations |3G | 

±dCdC ddCdC 

di'&f~~[h^'di~d9^'d9^' ^ ^ 



From ([3]) with ([5]) we obtain 



We now specialize the voltage V{t) to 

Vit) ^ Vdc - KcCos(r!t), (8) 

where Vdc > is the dc part of the voltage. Vac > is the ac part of the voltage, and ft is the angular frequency of 
the ac part of the voltage. We also define the dimensionless control parameters 

rj = ^ (9) 



and 



the unit of time 



al ' 



(10) 



*° = ^ (^1) 



and the unit of length 



3 



such that 



and 



T = — 

to 



(13) 



(14) 



are the dimensionless radius and time, respectively. In terms of p and t, with ^ - (jl4p . the dimensionless equivalents 
of ® and © are 



p+ [l-2^cos(2r)]^ = 0, 



+ I3[l - 2r/cos(2T)] 



where the dots now indicate differentiation with respect to r. 



sin(6i) 



0, 



(15) 
(16) 



III. LINEAR STABILITY ANALYSIS 



It is not obvious that the equations of motion (|15l) and (|16l) have stable solutions. Indeed, we will see in this and 
following sections, that both rj and /3 need to satisfy certain conditions in order for (|15p and to exhibit stable 
solutions. To find the condition for /?, we perform a linear stability analysis. We emphasize that expanding the 
equations of motion to linear order is not an uncontrolled approximation. In the absence of pathologies - as is the 
case here - it is an exact method for assessing the stability properties of the equilibrium solution of the equations of 
motion jsij . 

To first order in 9, we may replace cos{9) by 1 to obtain 



[1 - 2r;cos(2T)]- 



1 



0. 



(17) 



Thus, to linear order, the p equation decouples from the 6 equation. This fact substantially simplifies the stability 
analysis. However, the 9 equation remains coupled to the p equation since even in linear order the 9 equation contains 
p. To lowest order the Fourier expansion of the exact solution p{t) of p?)) is given by 



p(t) = Po - Pi cos(2r). 



(18) 



where po and pi are constants. Inserting into ((T7)) . replacing cos^(2r) terms by their average value 1/2, and 
neglecting {a + l)pi/pQ with respect to 277, the p equation (fT7)) is fulfilled if 



PO 

Pi 



a + 1\ °+ 



(19) 
(20) 



2°+2 {a + l)°+2 ?7=+2 

Inserting (|18l) into (jl6p . keeping only terms linear in pi, and replacing once more all cos^(2r) terms by 1/2, we obtain 

, /3 r/ 1 



Po 



a + l 



- 2r? 



1 - 



cos(2t) } sm{9) = 0. 



(21) 



Furthermore, since we are focusing on small oscillations, we may safely replace sin(0) by 9. In addition, as we will see 
below, rj is usually quite large (77 }t 10) so that we may neglect the term a/[2r]'^{a + 1)] with respect to 1. Together, 
this yields 



Po 



1 



a + l 



— 2i]cos{2t) 



= 0. 



This is a Mathieu equation [32|], whose canonical form is [34I 

X + [a — 2qcos{2T)]x = 0. 



(22) 



(23) 
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Comparing ([23]) with (|22)) we see that 



q 



Po' 



(24) 
(25) 



It is weU known [33 | that stable solutions of (I23p exist only in certain two-dimensional regions of the {q, a) parameter 
plane [S^l- Because of po > and the definitions of a and /3 in ([1]) and p^ . respectively, we have a > 0. In this case, 
and up to second order in q, the first, and most important, stability region of the Mathieu equation is bounded by 



< a < 1 



q 

q - ^■ 



With (EH) and ^ we have 



(a + 1)?7 

Using this in (|26p . we obtain the following quadratic inequality for q: 

1 

(a + 1)?7 

Keeping only terms up to first order in l/[(a + 1)?]], the solution is 



1 



< 0. 



where 



q < qsiv), 



qsiv) = (-4 + 2V6) - 



{a + 1)7] 



V6 
3 



= 0.9 - 



0.734 
(a + 1)77' 



With ([T9|) and ([25|) this yields the stability criterion 



13 < qsiv) 



Pa 



qsiv) 



a + 1\ °+ 



(26) 
(27) 

(28) 

(29) 

(30) 
(31) 

(32) 



In summary, a range of [3 values exists that leads to stable 9 oscillations. This result is significant. It means that 
as long as stable solutions of (fTT]) exist, stable solutions of the system (|15|) and (|16p can be constructed, which, in 
turn, implies stable confinement of dipoles. Therefore, before proceeding further, we need to investigate the stability 
properties of the radial equation (fTT)) . 



IV. PROPERTIES OF THE RADIAL EQUATION 

A pseudo-potential analysis [s^l is the most natural way to conduct a comprehensive stability analysis of the radial 
equation (jl7p . It is a powerful method that may be applied generally for the analysis of the motion of particles in 
rapidly oscillating force fields [l^l- According to this method, we start with ([T5)) . but instead of holding po fixed, we 
allow p{t) to execute slow oscillations around po- Therefore, we now write 

p{t) = R{t) ~ PiC0s(2t), (33) 

where, as in ()18p . pi is a constant. Equation ([33]) represents formally what actually happens physically: The exact 
motion p{t) is a superposition of a slow, large-amplitude motion R{t) and a fast, small-amplitude motion proportional 
to cos(2t) [30| . The slow, large-amplitude motion is known as the macro-motion; the fast, small-amplitude motion 
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is the micro-motion. Apparently, because it determines the cycle-averaged trajectory of a trapped particle, it is the 
behavior of R{t), which determines the stability properties of the trap. 

Inserting (pS)) into (fTTl) . keeping only terms up to first order in pi, and, as we did before, replacing cos^(2r) by its 
average value 1/2, we obtain 



R 



1 



1 



ri{a + l)pi 

Jla+2 



2?7- + } cos(2t) = 0. 



This equation is fulfilled if, separately, 



R 



R" 



1 



ri{a + l)pi 
Pi 



277 - (a -1-1) 



R 



= 0, 
= 0. 



Neglecting the term (a -I- l)pi/ R in (|36|) with respect to 2ry, equation (p6|) yields 

■n 



Pi 



2i?"+i ' 



Inserting this result into (|35|) yields 



R 



where 



V,s{R) 



dVesjR) 
OR 



1 



(34) 

(35) 
(36) 

(37) 

(38) 
(39) 



is known as the effective potential fso'] , the pseudo-potential [s^ , or the ponderomotive potential in atomic ^33] and 
plasma [sH] physics. Apparently, the motion of the point dipole, on average, behaves as if the dipole were under the 
influence of a force that is the gradient of the pseudo-potential Veg{R). This is a powerful result that allows us to 
determine the stability of the motion of the point dipole according to physical reasoning: If Vcg{R) has a potential 
minimum, the point dipole will execute stable, oscillatory motion inside of the potential minimum; the point dipole 
is trapped. If VcS does not have a potential minimum, the point dipole will escape to infinity, and no stable trapping 
is observed. A potential minimum occurs at points Rq for which V^g{RQ) = and V],'^(i?o) > 0. Indeed, V^ff{Ro) = 
has a single solution for 



Po, 



(40) 



where po is defined in (jl9l) . Apparently, linear stability analysis and pseudo-potential analysis are consistent in that 
they predict the same equilibrium radius po. Moreover, we have 



Kff(i?o) - 



2(a + 2) 
(a + 1)7?2 



> 0, 



(41) 



i.e., Rq — po, in fact, corresponds to a pseudo-potential minimum for all control parameters rj. At this point we have 
to remember that in our pseudo-potential analysis we neglected higher-order terms in pi/R. Thus, we expect the 
pseudo-potential analysis to be trustworthy only for pi/Ro 1, which implies 



77 > 



1 



1 



(42) 



Physically, the condition pi/Ro <C 1, or, equivalently, Rq ^ pi means that the amplitude pi of the fast oscillations 
around Rq needs to be small enough with respect to i?o such that the dipole does not hit the field-generating source 
at r = 0. This, however, would indeed occur for Rq < p ^ an d lead to the ejection of the point dipole from the 
trap. Indeed, there is "trouble" for small 77. It is known [3a] that for a — the radial equation (IT7|) exhibits a 
period-doubling bifurcation at rji w 3.12 and further period-doubling bifurcations, following the Feigenbaum scenario 
[31| . at rjj+i < rjj, j = 1,2,..., terminating in a transition to fully developed chaos at 7700 w 2.91. This means that 
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in rjj+i < T] < rjj, j = 1,2, .. ., the lowest Fourier component of p(r) is not cos(2t), as in (|18l) . but cos(2t/2-'). This 
is not taken into account in ([T51) . As a consequence, our pseudo-potential analysis is not valid for ij < iji. Detailed 
calculations show that (i) the period-doubling scenario persists for a > 0, and (ii) that ?yi(a) is a slowly, monotonically 
decreasing function of a, ranging from rjila — 0) « 3.12 to r]i(a — 10) ~ 2.41. It is important to emphasize that the 
occurrence of bifurcations does not preclude stable trapping. Indeed, stable trapping, in principle, may be achieved 
in the entire range [rjac < r] < iji]- Thus, the existence and location of the predicted bifurcations may be explored 
experimentally. However, since in this region, as mentioned above, the point dipole is already very close to the 
field-generating source(s), it may be very difficult to trap dipoles experimentally in the bifurcation regime. 

Another source of trouble for the pseudo-potential analysis are resonances, which occur for rj > rji. Resonances occur 
whenever the frequency of oscillations in the pseudo-potential well V^s is a rational multiple of the driving frequency 
of the trap. Luckily, only the lowest order resonances are "dangerous" . They occur whenever the low-frequency 
oscillations of R{t) in ([55]) are 1/3 or 1/2 of the driving frequency (= 2) in (|17p . Although the pseudo-potential 
analysis, incorrectly, predicts stability in these resonant cases, it is not entirely useless. Since rj > r]i, the conditions 
for the computation of the pseudo-potential are met. Therefore, there is nothing wrong with computing the frequency 
of slow oscillations of the dipole in the pseudo-potential VcS and, on the basis of this result, to predict the rj values at 
which the resonant instabilities occur. Since the frequency w of small oscillations around the equilibrium radius po is 



w - V^cff(Po), (43) 

the instability at = 2/3 is predicted to occur at 



l2{a- 


f 2) 


a-\ 


- 1 


'2(a- 


^2) 


a4 


-1 



= ^/^T^' (44) 
and the instability at a; = 1/2 is predicted to occur at 

V a -I- 1 

where we used the expression (|1T|) for V^g{po)- For a = 1, we predict r]^ = 3^/3/2 « 2.6 and rjl = 2^/3 « 3.5 in fair 
agreement with the numerically obtained values r]^ « 3.16 and rjl « 3.86, respectively (37| . More important than the 
numerical proximity of these special values is their behavior as a function of a. Equations (j44l) and (|45|) show that 
both 773 and rj^ are bounded and never exceed = 4. 

Summarizing the results of this and the previous section, we obtain the following general result: In the absence 
of gravity, stable trapping of point dipoles may be achieved for all ?/ > 5. Even for 77 < 5, ranges of i] values exist 
for which stable trapping is possible. The 77 < 5 regime exhibits rich dynamics, but is less accessible to analytical 
analysis. 



V. NUMERICAL INVESTIGATION OF STABILITY 



In this section we back up the analytical calculations with numerical simulations. We solved the system of equations 
(fTSl) and (fT6|) numerically for rj = 100, /3 ~ 0.59, and a — 2 for the following 5929 initial conditions: p(r = 0) = 
13+jx0.2, j = -5,...,5, p{t = 0) = fcxO.Ol, fc = -3, . . . , 3, 6'(r = 0) ^1x0.03, 1 ^ -5, . . . ,5, 9{t = 0) = mx 0.003, 
m = —3, . . . , 3. Since the solutions of the system of equations (|15l) and ([TS]) are well behaved, no special numerical 
integrator needs to be chosen. We chose a simple 4th order Runge-Kutta integrator [l^l, which yields completely 
converged, numerically exact solutions for any prescribed accuracy (we achieved a relative accuracy of better than 10~^ 
in all dynamical variables). Using the Runge-Kutta integrator, each of the 5929 initial conditions was propagated 
forward in time, and for each trajectory it was verified graphically that it stayed bounded over a time interval 
At = 5007r. 

A representative molecular trajectory is illustrated in Figs. [T] and [H The trajectory picked is the one with initial 
conditions p(0) — 13, p(0) = 0, 6{0) = 0.01, and 9{0) — 0. Figure [1] (a) shows the time-evolution of p(r) over the 
time interval < t < 5007r. During this time interval, the macro- motion performs about three and a half oscillations, 
ranging between pmin = 9.8 and Pmax — 13. Thus, the oscillations are bounded in the p direction. It is important 
to notice that the micro-motion is very small. In fact, the micro-motion amplitude is so small that it is not resolved 
on the scale of Fig. [T] (a) and manifests itself as a fine blur of the line. This is an important consistency check, 
since the pseudo-potential method is valid only if the micro-motion amplitude is small compared to the oscillations 
of the macro-motion. This is clearly born out in Fig.[T](a). Since rj is large, the small amplitude of the micro-motion 
observed in Fig.[T](a) is consistent with its analytical prediction ([^0)1 . 
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FIG. 1: Molecular trajectory calculation of (a) p(r) and (b) 9{t) for initial conditions p(0) = 13, p(0) — 0, ^(0) = 0.01, and 
8{0) = 0. The blurring of the line in (a) results from the micro-motion, which is not resolved on the scale of (a), and, as (a) 
shows, is very small. Panel (b) shows bounded oscillations of S as a function of time. Together, (a) and (b) illustrate the 
stability of the trap. The calculations are numerically exact (converged) solutions of the equations of motion (|15|l and (|16|) . 
No approximations were performed. 

Figure [2] illustrates the stability of the representative molecular trajectory in phase space. This is a more stringent 
test of stability than the p(t) and 0{t) trajectories displayed in Fig. [1] since the trajectory, e.g., might be stable in 
p, but unstable in p. We show the molecular trajectory over the same time interval chosen in Fig. [1] i.e. At = SOOtt. 
Since the molecular trajectory is embedded in a 4-dimensional phase space (p(r), p(t), 9{t), 0{t)), we use the technique 
of Poincare sections [3l| to display the trajectory via six panels of two-dimensional phase-space projections. We choose 
the projections on the six possible two-dimensional coordinate planes (a) p vs. p, (b) 9 vs. p, (c) 9 vs. p, (d) 9 vs. p, 
(e) 9 vs. p, (f) 9 vs. 9. While the micro-motion is essential for generating the trapping forces, it is the macro-motion 
that determines the stability of the trap (see Sec. lIV[) . Therefore, suppressing the irrelevant micro- motion and focusing 
on the relevant macro- motion, we strobe the motion, i.e., we display the trajectory at the end of each cycle of the 
trap field at the discrete times tj = j x ir, j = 1, . . . , 500 (plot symbols in Fig. [5]). Figure [5] (a) shows clearly that 
the trajectory executes a stable oscillation in p direction. Figure [2] (b) shows that the motion is also bounded in the 
p-9 plane. Figure [2] (c) shows new information; we see that the motion is also bounded in 9 direction. Figures [2] (d) 
and (e) contain the new information that the motion is also bounded in p. Figure O (f ) shows that the trajectory 
executes stable oscillations in the 6 direction. The six panels taken together imply that there is no "escape direction" 
in phase space, i.e. the motion is globally stable. Plots akin to Figs. [T] and [5] were produced for each of the 5929 
initial conditions and overlaid one of top of the other. The resulting composite plots look qualitatively like Figs. [1] 
and [21 showing no data points outside a confined area. This way it is possible to summarily assess stability of all 
5929 molecular trajectories, without having to inspect them one-by-one individually. The net result is a numerical 
confirmation of stability of the trap in an entire four-dimensional phase-space volume. 

Comparing Fig. [T](b) with Figs. [2] (b), (d), and (f), one might perceive a problem: The 9 amplitude in Fig. [T](b) is 
larger than the 9 spread in Figs. [5] (b), (d), and (f). The reason for this, however, is straightforward. It was checked 
explicitly that the reduction of spread in Figs. [2] (b), (d), and (f) is due to the strobing of the motion at multiples of 
TT. Strobing does not display the motion over a continuous time interval, but samples the motion at discrete times. 
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FIG. 2: Six stroboscopic sections of the 4-dimensional phase space of the molecular-dynamics trajectory of Fig.[T]for 500 cycles 
of the trap field. Each point in the figure represents the respective values of the trajectory at the end of each cycle, (a) p vs. 
p, (b) 9 vs. p, (c) 6 vs. p, (d) 9 vs. p, (e) 9 vs. p, (f) 9 vs. 9. The six sections show clearly that (i) the molecular trajectory is 
bounded in all four phase-space dimensions, and (ii) executes bounded oscillations in both p and 9 directions. 



This explains the small reduction in spread. 

The numerical check performed amounts to a stability survey of a four-dimensional phase-space volume surrounding 
the point {p{t = 0),p(r = O),0(t = 0),^(r = 0)) with a resolution of 5929 initial conditions that all lead to stable 
solutions of the system of equations (|15|) and ()16|) . The composite phase-space plot produced by overlaying all 
5929 individual phase-space plots reveals stable radial oscillations in the p interval 9 < p < 19 and the 9 interval 
-0.3 <e < 0.3. 

Since dipole forces are usually much smaller than corresponding forces on isolated electric charges, we anticipate 
that consideration and inclusion of gravitational forces may be important for a proper description of the dynamics of 
trapped dipoles according to the new trapping scheme. We will address the influence of gravity next. 



VI. BEATING GRAVITY 



In previous sections we established that, in the absence of gravity, control parameter combinations (jy,/?) exist that 
result in the stable, permanent confinement of point dipoles. Since in this paper we are focusing on electric fields 
of the form ([l}, which vanish for large r, it is obvious that in certain trapping regimes gravity may overwhelm the 
confining effect of the electric field and result in a destabilization of the trap. Therefore, in this section, we study the 
effects of gravity and state the conditions that result in stable trapping in the presence of gravity. 

In SI units, and according to (p9)) . the pseudo-potential Vcs{r) is given by 



1 



where p = r/lo is defined in (|13p . In the presence of gravity, Vcf[{r) is modified to Vcs{r) — )■ Vcs{r) — Mgr, i.e. 



tl 



(46) 



(47) 
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where 



with 



1 7f _ 

94 



7 



^0 



For stable trapping to occur, Vcff{p) needs to exhibit a potential minimum, i.e. we require 

vM = /(p) - 7 = 0, 

where 

1 + 



fip) 



P' 



A+i 2p'^°'+^ 



(48) 



(49) 



(50) 



(51) 



Since 7 > 0, equation (1501) has a solution only if the maximum of f{p) exceeds 7. From f {p) = we determine that 
the maximum of f{p) occurs at 



/2a + 3 



which results in the minimum condition 



2o + 3 

a + 2 / 2 \ 



7 < fiPmHK) -— 2C. + 2 1 D I o 

27^t;+2- V 2a + 3 



(52) 



(53) 



VII. DESIGN ALGORITHM 



At this point we know that the new mechanism is capable of trapping dipoles, but several conditions have to be 
fulfilled simultaneously. The following algorithm provides a systematic way for designing trap parameters that satisfy 
the various conditions. 



1. We start with the stability criterion ([32l) and use the definition of (3 in ([T0| to obtain the condition 



lo < 



qs{ri)aIp'S 
Mr] 



1/2 



(54) 



Denote by D the desired location of the radial equilibrium distance of the trapped dipole. The size of D may 
depend on many factors, prime among them the spatial extent of the field-generating electrodes. If, e.g., the 
field is generated by a wire, we obviously need to choose D > Twirc, where r^h-c is the radius of the wire. With 
D chosen, and 



D = Polo, 



we now have 



Pi 



/2 
1)2 



< 



El 
V 



With dig]) this implies 



f] > 



qs{ri) a{a + 1)1 



(55) 



(56) 



(57) 



Since the right-hand side of ([57)) . via qs{i]), depends on 77, we solve ([57)1 iteratively. Starting with qsif]) = 0.9, 
an 7] value satisfying ([57)) is usually obtained after only a few iterations of ([57]) . 
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2. With rj chosen, we may now use ([55]) to compute Iq according to 

, _ D _ D_ 



Since r] satisfies ([57)) . Iq, computed according to ([55]) . automatically satisfies ([M)) . 



(58) 



3. With r] chosen and Iq computed, we may now derive a condition for the drive frequency fl. From (|53p with (j49p 
and ini), we obtain 



n > 



8g 



(a + 2)/o 



2q + 3 

2a + 3\ 



1/2 



(59) 



4. With rj, Iq, and $7 known, we now use (|12p to determine Vdc according to 



27Q + 2 







Aapv 



(60) 



5. With all design parameters chosen, we may now perform an important verification step, the numerical solution 
of the system of equations 



p+ [l-2r;cos(2T)]^ -7 = 0, 



+ /3[1 - 2rycos(2r) 



sin(6') 



0, 



(61) 
(62) 



i.e. the system ([TS|) . (ITC)) including the effect of gravity in the radial equation. Only if the system (pT|) . (p^ 
exhibits stable solutions did we indeed construct a parameter set that leads to the stable trapping of point 
dipoles p in the field E. Failure of this verification step indicates that the trap parameters were chosen too close 
to the analytical estimates of the stability borders. Allowing for larger safety margins will cure the problem. 



VIII. FERROELECTRIC RODS 



Ferroelectric rods of mass M, trapped in the electric field of a metallic, spherical electrode of radius rc, provide an 
illustration of the new trapping mechanism that may be realized experimentally. To be specific, we study the case of 
thin rods of barium titanate, a classic ferroelectric material [ssl ] with a spontaneous polarization of Pg — 0.15 Cm~^ 
[39} and a density of 6g/cm'^ [i^]. The moment of inertia of thin, cylindrical rods is [41] 

I = ^Af6^ (63) 

where b is the length of the rods. We specify 6=1 mm, D = 2 mm, and z/ = Tc = 1 mm. According to step 1 of the 
design algorithm and qs = 0.9, we obtain ij > 17.8. We choose ij = 25, which satisfies (|57)) . According to step 2 we 
obtain Iq = 0.36mm. According to step 3 we obtain ft > 7811 s~^; we choose f2 — 1.5 x 10^ s^^. Step 4, then, yields 
Vdc = 18.9 V. The checks on ^ and 7 yield (3 = 0.78 < QspI/v = 1-1 and 7 = 4.8 x lO"'^ < 2{2/7y/^ /rj^^^ = 1.8 x lO'^. 
According to (j20p the micro- motion amplitude of the rod is pilo — 26.6 /im, which allows enough clearance for the 
rod to oscillate around r — D. The numerical check of stability according to step 5 was performed for rj — 2b, 
(3 = 0.78, a = 2, and 7 = 4.8 x 10""* over a time interval of At = 5007r for the following — 2401 initial conditions: 
(p(r = 0) = 5.3 + j X 0.03, p(r = 0) = fc x 0.003, 6I(t = 0) = / x 0.01, ^(r = 0) = m x 0.01), j, fc, /, m = -3, . . . , 3. It 
was verified graphically that all 2401 trajectories are stable. 



IX. DISCUSSION 



The new ac/dc trapping mechanism proposed here is fundamentally different from the saddle-point mechanism 
employed for the successful demonstration of electrodynamic trapping of polar molecules j23l . [2^ . The trap proposed 
here makes use of monopole fields, i.e. particles are trapped in close proximity to sources of the electric field. In 
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addition, while the pseudo-potential of higher-multipole traps generates a focusing force, the pseudo-potential of the 
trap proposed here generates a defocusing force, which attempts to drive the particle toward r — oo. Thus, the dc 
voltage is an essential ingredient for the trap discussed in this paper: the dc voltage generates the attractive force 
counterbalancing the defocusing effect of the ac electric field. Thus, the trap discussed here is akin to the dynamically 
stabilized inverted pendulum (42„ .43,] , a classic lecture demonstration in introductory mechanics courses. 

At a first glance, the new trap introduced here resembles a trap investigated several years ago by Riis and Barnett 



44] . However, there are several important differences, (a) The Riis/Barnett trap aims to trap molecules in the 



ground state [44|, whereas the new trap employs a dynamical trapping mechanism in the radial and angular degrees 
of freedom. Quantum mechanically, trapping in the angular degree of freedom corresponds to a wave function in 
which the trapped particle finds itself in an angularly localized wave packet that does not correspond to the rotational 
ground state of the trapped particle, (b) For the Riis/Barnett trap gravity is essential [3, EH, whereas the new trap 
works even without gravity (see Sec. [V]). (c) In contrast to the Riis/Barnett trap, which does not use a permanent 
dipole moment, a permanent dipole moment is essential for the new trap to work, (d) No stable trajectories have 
been found [i^ for the Riis/Barnett trap. Far from being trivial, the reason why the Riis/Barnett trap does not work 
is interesting and instructive: The micro-motion in the Riis/Barnett trap is too large (45*. Contrary to the dynamics 
exhibited by the Riis/Barnett trap, the micro-motion in the new trap is very small (see Fig. [Ij, explaining why the 
new trap works. In our case, the small micro-motion (see Fig. [T]) is the reason why the pseudo-potential method is 
valid and leads to stable trapping. The importance of a small micro-motion for stable trapping and for the validity 
of the pseudo-potential method was also pointed out by Cornell [1^ . 

Our analysis did not include the possibility that the trapped dipoles may have orbital angular momentum, i.e., the 
possibility of revolving around the field-generating electrodes. Since orbital an gula r momentum is not an essential 
part of the trapping mechanism, like it is, e.g., in the classic Kingdon trap p4l[l3.ll5j. this is not an essential omission. 
On the contrary: It is one of the strengths of the new trap that it works for zero orbital angular momentum. In 
fact, in the presence of a cooling mechanism, zero orbital angular momentum will be achieved automaticall y. Z ero 
orbital angular momentum is a problem for traps that need orbital angular momentum to work (see, e.g., (20l l2l|): it 
is not a problem for the new trap. It has been checked that the new trap does work in the presence of orbital angular 
momentum. Some aspects of this result are not surprising. In the case a — 2, e.g., orbital angular momentum creates 
a centrifugal force with the same radial dependence as the dc electric force. Thus the effects of angular momentum do 
not structurally change the equations of motion and can be counterbalanced by an adjustment of ac and dc voltages. 
The fact that the orientation of the dipole on its trip around the electrode adiabatically follows a near- alignment 
with small oscillations around the orbital radius vector is more of a surprise. This mechanism, however, is not new; 
it forms the dynamical basis of the Sekatskii/Schmiedmayer trap 21 L which has already been demonstrated to 
work experimentally (47| . In fact, the Sekatskii/Schmiedmayer trap 20r 21 1 is a limiting case of the new trap in the 
presence of orbital angular momentum but zero ac electric field. Since the stabilizing repulsive force counteracting 
the attractive dc field is provided by orbital angular momentum for zero ac field, and by the repulsive pseudo-force for 
zero orbital angular momentum, a "chirping mechanism" may be envisioned: As the particle loses angular momentum 
due to cooling or damping, an ac voltage may be switched on, providing the repulsive force that keeps the trapped 
particle(s) from crashing into the electrode(s). The final state is then the zero orbital angular momentum state, fully 
stabilized by the ac field according to the new mechanism. 

The trapping mechanism discussed here is also different from the mechanism of optical tweezers (48l - [50j . Optical 
tweezers work astonishingly well for small particles ranging from atoms and molecules to biological cells and even 
buoyant, small organisms ^1,|4^. However, in contrast to the scheme proposed here, optical tweezers may not be 
able to hold a mm sized rod of ferroelectric against gravity, and the frequency of optical tweezers (~ 10"'^^ Hz) is orders 
of magnitude larger than is required according to the new trapping scheme. Thus, differing in mechanism and range 
of application, each of these two mechanisms has its own optimal niche of applicability. 

In our analysis above, we used E l/r" and examples were presented for a = 2. Inverse powers is a natural choice 
(wires, e.g., generate a — 1- spheres generate a = 2; oppositely charged spheres, e.g., generate a = 3 for large r), but 
not necessary. Superpositions of such fields as well as many other field configurations may lead to stable trapping, 
provided the corresponding stability criteria are worked out and are satisfied. 

For the time-dependent part of the electric field of the trap we chose a cos(f2i) drive (harmonic drive). This, again, 
is not essential. Any drive envelope, such as, e.g., rectangular or impulsive ("5 kicks"), can be accommodated and 
may yield stable, trapped solutions following proper adjustment of the /3, 7, and rj criteria. Impulsive drives, in 
particular, offer the possibility of representing the dynamics of the dipole in the form of iterated mappings 31]. This 
considerably simplifies the analytical and numerical analysis of the trap. 

Concerning the case of trapped barium titanate rods, we have to be careful that the applied electric fields do not 
exceed the coercive field strength, i.e. a field strength that, applied in the opposite direction, would result in a reversal 
of the direction of the permanent polarization of the rods. Since, according to V"ac = 2riVdc, and Vdc = 18.9 V in 
our example, the ac electric field at the position D = 2mm is i?ac = rcVac/D^ = 2400 V/cm. Since typical coercive 
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fields for barium titanate are on the order of lOkV/cm [5l|, the field in our example is tolerable. However, if actual 
experiments are contemplated, it may be advisable to optimize the trapping parameters and/or the nature of the 
trapped particles. Rods made of croconic acid crystals (S^l, e.g., may be easier to trap than the barium titanate rods 
discussed in Sec. I Villi 

Another concern is the large dielectric constant of barium titanate, which may result in an appreciable induced 
dipole moment. At room temperature, e^, the relative dielectric constant, may reach values of tr ~ 5000 [s^. In this 
case the induced polarization is Pi„d = ^o^rEac = 1-1 x 10~^C/m^ <C Ps = 0.15 C/m^, where eo is the permittivity 
of the vacuum. We conclude that induced polarization is not an important effect. Even if for some trapped particle 
species this effect should be large, the induced polarization may be fully compensated for by adjusting Vac and Vdc- 

While this paper focuses on the case of a single stored dipole, we may also study the case of many, simultaneously 
stored dipoles. In the presence of a cooling mechanism, such as laser cooling [5j] (for atomic-sized particles) or 
buffer-gas cooling [H?] (for atomic-, nano-, and macroscopic-sized particles), trapped particle densities may become 
large enough that dipole-dipole interactions are important. In the case of molecules this may allow us to study new 
ferroelectric phases of dilute neutral (trapped) molecular gases. In the case of small ferroelectric particles, in analogy 
to the geometric arrangements of trapped Coulomb crystals [T6l-[l8^. a geometrically ordered dipolar phase may be 
generated and studied. This may result in an interesting contribution to the emerging field of granular materials [5^ . 



X. SUMMARY AND CONCLUSIONS 



This paper presents a novel electrodynamic trap for the stable confinement of neutral particles with permanent 
electric dipole moments. The mechanism of the trap differs profoundly from currently demonstrated or proposed 
neutral-particle traps. While conventional electric traps are based on energy shifts in selected quantum states and 
employ higher-multipole fields as an essential ingredient of trapping, the new trap uses a superposition of monopole 
fields, i.e. an attractive, static dc electric field and an ac field, whose net effect results in strong repulsion. If certain 
conditions on the strengths of the dc and ac fields are met, stable trapping results at the point where attractive and 
repulsive forces are in equilibrium. A mechanical analogue of this trapping mechanism is the dynamically stabilized 
inverted pendulum. Stability of the trap is proved with the help of an analytical stability analysis, supplemented with 
detailed numerical molecular trajectory simulations. The effect of gravity is considered and included. The example 
of thin barium titanate rods indicates that the trap works in practice. It is possible that the trap also works for polar 
(macro) molecules. Although there are no magnetic charges, and magnetic fields are divergence-free, it is possible 
that close to a magnetic pole as a field generating "source" a superposition of static and oscillating magnetic fields 
may result in a magnetic analog of the trapping mechanism discussed here. 
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